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Abstract: 

In this work we provide all the missing ingredients for calculating parity-violating 3-jet observables in 
e^e^ collisions at next-to-leading order (NLO) including the full quark mass dependence. In partic- 
ular we give explicit results for the one-loop corrections and for the singular contributions from real 
emission processes. Our formulae allow for the computation of the NLO 3-jet contributions to the 
forward-backward asymmetry for heavy flavours, which is - and will be - an important observable for 
electroweak physics at present and future e+e^ colliders. 
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I Introduction 



The measurement of the forward-backward asymmetry Apg for ^-quark production in e^e^ coUisions 
at the Z peak provides one of the most precise determinations of the weak mixing angle sin"6eff. It can 
be determined with an error at the per mille level if the computations of Apg are at least as accurate as 
the present experimental precision of about two percent This implies, in particular, that perturbative 
and non-perturbative QCD corrections to the lowest order formula for Apg must be controlled to this 
level of accuracy. 

As far as perturbative QCD corrections to this observable are concerned, the present knowledge is as 
follows. The order contributions had been computed first for massless and then for massive quarks 
in ref. [Q] and ref. ||3|,^, respectively. These calculations, which used the quark direction for defining 
the asymmetry were later modified [§-0] to encompass the asymmetry with respect to the thrust axis, 
which is more relevant experimentally. The coefficient of the order correction was computed first 
in ref. for massless quarks using the quark axis definition, by numerical phase space integration. 



This result was recently corrected by a completely analytical []9|| and by a numerical calculation [|10p. 



In the latter reference the second-order corrections were also determined for the thrust axis definition. 

For future applications to measurements at proposed linear e^e^ colliders the attained theoretical preci- 
sion so far will not suffice. The forward-backward asymmetry will be a key observable for determining 
the electroweak couplings of the top quark above threshold - a fundamental task at a linear collider as 
far as top quark physics is concerned. Needless to say all radiative corrections must be determined for 
non-vanishing top quark mass. Moreover, there may be the option to run such a collider at the Z peak, 
aiming to increase the precision attained at LEPl/SLC significantly. For instance one may be able to 
measure Afb for b and c quarks to an accuracy of order 0.1 %. (For a recent discussion, see ref. Jill].) 



Here we are concerned with the second-order QCD corrections to Afb for massive quarks Q. It is 
convenient to obtain this asymmetry by computing the individual contributions of the parton jets [||]. 
This break-up may also be useful from an experimental point of view, as the various contributions can 
be separately measured [12-14]. To order the asymmetry is determined by the 2-, 3-, and 4-jet final 
states involving Q. The 4-jet contribution to Afb involves the known parity- violating piece of the Bom 



matrix elements for e^e^ 4 partons involving at least one QQ pair (see, e.g. ref. Q15|]). The phase- 
space integration of these matrix elements in order to get (AFB)4-jet for a given jet-resolution algorithm 
can be done numerically in a rather straightforward fashion. The computation of (AFB)3-jet is highly 
non-trivial. Within the phase space slicing method [[1^18], which we will use here, it amounts to 
determining first the parity- violating piece in the fully differential cross section for the reaction e^e^ — > 
3 resolved partons to order a^. This will be done in this paper. (For the definition of 'resolved partons' 
see section IV.) In addition the contribution to (AFB)3-jet from e^e^ —>■ 4 resolved partons is needed. It 
can be obtained from the known leading order matrix elements by a numerical integration in 4 space- 
time dimensions which incorporates a jet defining algorithm. This will be discussed elesewhere. The 
computation of (AFB)2-jet for heavy quarks is beyond the scope of this work and remains to be done in 
the future. 

Our paper is organized as follows. In section we review the kinematics of the reaction and leading- 
order results, and define the 3-jet forward-backward asymmetry. In section [11 we discuss the calcu- 



lation of the virtual corrections. In particular, the ultraviolet divergences, the infrared divergences, 



and the renormalization procedure is discussed. In section |V| the singular contributions from the real 



corrections are calculated and the cancellation of the infrared divergences is checked. In section ^ we 



1 



summarize our results and end with some conclusions. 



II Kinematics and leading-order results 

Tlie leading order differential cross section for 

e+{p^)e-{p^) ^ Q{kQ)Q{kQ)g{kg) (II.l) 
can be written in the following form^: 

„ , , , = (l+cos2i3)Fi + (l-3cos2i^)F2 + cosi3F3 

d<^dcos-&dxdx 4 (47i)3 J ^ \ ) i J 

1 

+ sin2'dcos(j)F4 + sin i^cos2(j)F5 + sin'dcos(j)_F6 > , (11.2) 



with 



471 fy 

C5p, = a[e^e- ^^^ix-) = (II.3) 



In eq. (11.2), denotes the angle between the direction p_ of the incoming electron and the direction 
kg of the heavy quark Q. The angle (j) is the oriented angle between the plane defined by and Q 
and the plane defined by the quark anti-quark pair. The functions Fi depend only on the scaled cm. 
energies x = 2k ■ kg/s and x = 2k ■ kg/s, and on the scaled mass square z = m^/s of the heavy quark 
and anti-quark, where k = + and s = k^. Note that the leading order differential cross section 
for a final state with n>4 partons can be written in an analogous way in terms of the two angles "d, (j) 
and 3n — 7 variables that involve only the final state momenta. In higher orders absorptive parts give 



rise to three additional functions F7 g 9 (cf. ref. []19|]). The four functions Fi 2,4,5 are parity even. In 



particular, the 3-jet production rate is determined by Fi. In the case of massive quarks next-to-leading 
order results for Fi were given in references [2^27]. The two functions F3 and F^ are induced by 



the interference of a vector and an axial-vector current. In particular, using the quark-axis, the 3-jet 
forward-backward asymmetry is related to F3 in leading order in the following way: 



('^Fb) 3-jet 



/o ^ ^I''nd<\> //)(,■„,) dxdx^^^^^^- - /"i J COS 1^ /^^^ fifcj) ^^^■^ ^c^rfcofwf 
/- 1 COS 1^ /^'^ dxdx^j^j^^^- 

3 JD{y,,,)dxdm{x,x) 



^ SDiy,,Adxdm{x,xy 



(11.4) 



where D(ycut) defines, for a given jet finding algorithm and jet resolution parameter y^ai, a region in the 
{x,x) plane. Analogously, the function induces an azimuthal asymmetry 

_ d<^ J\d cos^ /^(y^^^) d^d^ d^dcosldxdx - In/2 d(^!\dcos-& /^(y^,^^) dxdx ^^^J^%^^^- 

r'^/^^Afi ^^r^cA r ^^^^ da I r't ^,.,^cAr ^^^^ da 



lo d<\>J_yd cos 1^ /^(y^^^) dxdx ^^^JX,^^. + /;/2 d<^J_^d cos i3 /^j,,^.^,) dxdx^^^j^^^^^ 
3 lD(y,,,)dxdxF6ix,x) 
^ lDiy,,,)dxdxFiix,x) 



(11.5) 



We neglect the lepton masses and do not consider transversely polarized beams. 
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Note that final states with two partons do not contribute to A,^. The electroweak couplings appearing in 
F3 6 can be factored out as follows: 



(11.6) 



Ineq. (II.6), 



J - rf 



5a ''3 ' 

X 



-2e/sin''ew, 

1 



4 sin 9iv cos^ Q^/ s — m2 + imzTz 

= (11.7) 

where is the third component of the weak isospin of the fermion /, Qw is the weak mixing angle, 
and Xip denotes the longitudinal polarization of the electron (positron). In next-to-leading order in a^, 
additional contributions to F3 6 with electroweak couplings different from those in eq. ( |II.6| ) are induced 
which we will not consider here. They are either proportional to Imx and thus suppressed formally in 
the electroweak coupling or generated by the triangle fermion loop diagrams depicted in figure 1(1), 
l(m). The contribution of the triangle fermion loop is gauge independent and UV and IR finite. It was 
calculated some time ago in ref. 



The functions F3 6(x,x) may be expressed in terms of functions /zg, hj which appear in the decomposi- 
tion of the so-called hadronic tensor as performed for example in references [p^31|. 



^3 = ^ -4z/j6(x,x) + ^/ - 4zcos-&QQhT{x,x) 
F(, = -^\/^^^sini^gg/j7(x,f), (II.8) 

where -^qq is the angle between Q and Q in the cm. frame and we have 

l{\-x-x + lz)+xx 
cos-d^Q = (II.9) 
Vx2-4zVx2-4z 

It would seem pointless to trade ^3^5 for /jg 7 if it were not for the relation (which follows from CP 
invariance): 

hT{x,x) = —h(,{x,x). (11.10) 

The objective of this paper is to provide explicit expressions for the virtual corrections and for the 
singular parts of the real corrections to the function h^. These corrections make up the complete 
contribution from three resolved partons at NLO, which is ultraviolet and infrared finite. We 

regulate both infrared und ultraviolet singularities by continuation to J = 4 — 28 space-time dimensions. 



In this context one may wonder whether the above kinematics relations, in particular eq. ( |II.8| ) have to 
be modified for d y^4. This is however not necessary if we consistently keep the momenta of electron 
and positron as well as the polarization vector of the photon/Z-boson in J = 4 dimensions throughout 
the calculation. Note that this prescription still corresponds to the usual dimensional regularisation. 
In addition this prescription is compatible with the 't Hooft-Veltman prescription for 75 as we discuss 
later. 



3 



We finish this section by giving the Born result for h(,in d = 4 — 2e dimensions: 



/jg = 167ia,(A^'^- l)B[ 2x- {xx + r + 2-4x)e-4 



\-x 



with 



B 



{\-x)(\-xr 



and the scaled gluon energy 



(11.11) 



(11.12) 



(11.13) 



Note that the terms proportional to £ in eq. ( [I.ll| ) depend on the prescription used to treat 75 in d 
dimensions. To derive the above equation we have used the prescription 7^75 ^E^pyS V^Vy^ [ p2| , 
33]. We will discuss this issue in more detail in the next section together with the ultraviolet (UV) 
renormalization. 



Ill Virtual corrections 








(e) 





(1) 



(m) 



Figure 1 : One-loop diagrams to Z ^ QQg- 



In this section we will give analytic expressions for the virtual corrections to /jg, isolate the ultraviolet 
and infrared singularities, carry out the renormalization procedure, and discuss the treatment of 75. 
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We work in renormalized perturbation theory, which means that the bare quantities (fields and cou- 
plings) are expressed in terms of renormalized quantities. By this procedure one obtains two contri- 
butions: one is the original Lagrangian but now in terms of the renormalized quantities, the second 
contribution are the so-called counterterms: 

i:(*Po,Ao,mo,go) = L{^R,AR,mR,gR) + Lcti^R,AR,mR,gR). (III.l) 

The first contribution yields the same Feynman rules as the bare Lagrangian but with the bare quan- 
tities replaced by renormalized ones. In the following we renormalize the quark field and the quark 
mass in the on-shell scheme, but treat the gluon field and the strong coupling in the modified minimal 



subtraction (MS) scheme [34, 35], taking into account the Lehmann-Symanzik-Zimmermann residue 
in the S-matrix element. Note that the conversion of the on-shell mass to the frequently used MS mass 
can be performed at the end of the calculation. The quark mass in the on-shell scheme will be denoted 
in the following simply by m for notational convenience. 

While in the MS scheme the renormalization constants contain only UV singularities, in the on-shell 
scheme they contain also infrared (IR) divergences^. Although at the very end all the divergences 
cancel it is worthwhile to distinguish between UV and IR singularities so that one can check the UV 
renormalization and the IR finiteness independently. 

To start with let us discuss the contribution of the one-loop diagrams before renormalization, that is 
the contribution from L{^'R,AR,mR,gR). Although tedious the calculation of the one-loop diagrams 
shown in figure [l| is in principle straightforward. We do not discuss here the contribution from dia- 
grams 1) and m) as mentioned earlier. Without going into the details we just review some techniques 
used in the calculation of the remaining diagrams. To calculate the one-loop corrections we decided 



to work in the background field gauge [36,37] with the gauge parameter set to one. In this gauge the 
three-gluon vertex is simplified which leads to a reduction of the number of terms encountered in in- 
termediate steps of the calculation. To reduce the one-loop tensor integrals to scalar one-loop integrals 
we used the Passarino-Veltman reduction procedure [|38|]. From the diagrams shown in figure |l| one 
can immediately read off the scalar box-integrals one encounters in the calculation. These integrals are 

D^^' = J- [d^l (^^^)'''^' (III 2) 



° in^J (f + ie)((l + ko?-m^ + ie)((l + ko,? 



+ /8)((/ + kgY - m2 + 18) ((/ + kgg)^ - nfi + - kQY -m^ + ie) 



j~.d,sl.2 Y^d.sl.l 

^0 ~ ^0 



(III.3) 



(111.4) 



with kQg = kQ + kg. Note that only the real parts of the one-loop integrals contribute to /jg since we 
neglect terms in the electroweak couplings ~ Im^. For simplicity we will leave out in the following 
the prescript Re in front of the integrals. All the lower-point integrals are defined with respect to the 



box-integrals. We adopted the notation from Passarino and Veltman [ |38| ] in which capital letters A, 
B,. .. are used to define one-point, two-point,. . . integrals. The kinematics is specified by denoting 
as argument the propagators which are kept with respect to the defining box-integrals. For example 
Co(l,2,3) which originates from Dq is given by 

<^o( 1 ' 2, 3) = jdH ^ .^^^^^ ^ ^^^^ ^^^^ + + ;e,)2 - m2 + /£) • 

^ Note that the cross section at hand is inclusive enough so that no uncancelled collinear singularities survive. For 
simplicity we thus call soft and collinear singularities just IR singularities. 
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We have substituted the scalar box integrals in d dimensions by box integrals in <i + 2 dimensions using 
the relation 

Di = -^D^^\ (111.6) 
where is the coefficient of g^v in the decomposition of the four-point tensor integral Djjy (cf. eq. 



(F.3) in ref. []38|]), which in turn can be expressed as a linear combination of the scalar box integral 
Dq and scalar triangle integrals Co in d dimensions. Since the box integrals are finite in 6 dimensions, 
all IR singularities reside in triangle integrals after this substitution. As a gain, the coefficients of the 
triangle integrals become much simpler. This is clear, since the coefficient that multiplies the sum of 
all IR singular terms has to be proportional to the Born result. 

We may group the diagrams shown in figure |l| according to their colour structure in terms of the 
colour matrices Ta of the SU(N) gauge group. In particular, we distinguish between contributions 
which are dominant in the number of colours ('leading-colour') and contributions which are suppressed 
in the number of colours ('subleading-colour'). The different colour structures are separately gauge 
independent. The diagrams a, b, c are proportional to N{Ta)cQCQ where a is the colour index of the gluon 
andcg (cg) the colour index of the (anti-) quark. The diagrams h,j, k, 1 are proportional to l/N{Ta)cQCQ- 
They are purely sub-leading in the number of colours. The diagrams d, e, f, g are proportional to 
(A'^ — 1 /N){Ta)cQCQ- These diagrams contribute to the leading as well as to the subleading parts. 

To organize the next-to-leading order virtual corrections /j^"^ ' '^^^ for h(, it is convenient to separate 
them into different contributions. First we split the result into a contribution from the loop diagrams 
(figure |l|) and a contribution from the renormalization procedure K^^ ' . The contribution from the 
loop-diagrams is further decomposed into a UV divergent {K^^ ' '^" ), a IR divergent {K^^ ' '^^^') and 
a finite part (/jg'" ' *'" ). We thus arrive at the following decomposition: 

; virt., NLO 7 virt., UV div. , > virt., IR div. , ; virt., fin. , > virt., ren. /ttt ^^ 

K =K +K +K +K ■ (ni-7) 

We define the singular contribution by taking only the pole-part of the loop-integrals but keeping the 
8-dependence in the pref actors. 



III.l UV singularities 

All the UV singularities in a one-loop calculation appear in the scalar one- and two-point integrals. 
Defining the finite contributions A(m^), Bo{i,j) of these integrals by 

— ^ m - 
4TlfJ^ J £ 

m \ 1 ^ 



^o{iJ) = rr{^] -+Bo{i,j), (III.8) 



where rp is the usual one-loop factor: 



r(i+£)r2(i-£) 
r(i-2s) • 



the contribution from the pole-part of these integrals is given by 
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with 



and 



+ 8a?(A^2-l)C/7B8^1+;c-2x+(l-^)^ 



+ z 



lB{{f- - 2x){\-x) + (2 - 3x)(l -x)) +^1 



x{\—x) 

+ \6<xICfBz{—{xx -Ax + 2 + x^)- 2zXgB{x -x)) + 0{z^)\ 

^8 J 



+ 4zB(2. 



\-x 



■+Xe 



(III. 10) 



2{x^ — 3x — x'^x + 3xx + 2x^ — 4x + 2) —Xg{x'^ +xx — 4x + 2)e 



-2 



l-x 



{-5x-3x + 2x^ + 4 + x^ +xx)z^ , 



(1 —x)(x — 2x) , , 



(III. 11) 



(III. 12) 



The finite contributions involving A, Bq are given in section III.3. The cancellation of the UV singular- 
ities by the renormalization procedure will be discussed in sec, pl.4 . 



III.2 Soft and coUinear singularities 

To start with let us consider the soft and collinear singularities appearing in the leading-colour con- 
tribution. As mentioned earlier we use the box-integral in 6 — 2s dimensions rather than in 4 — 2e 
dimensions. The one- and two-point integrals encountered in this calculation and the box integrals in 
6 — 2s dimension do not contain IR singularities. As a consequence all IR singularities appear in the 
triangle integrals. Using the results for the integrals given in ref. |21] we see that only Cq(1,2,3) and 
Co(l,2,4) are IR divergent. The contribution from the two integrals is given by 



viit., IR div. 



lead, colour 



Using ref. [g!]] 



-gA^/j^°{(l-x)5C/,(l,2,3) + (l-x)5C^(l,2,4)}. (III.13) 



x)54(l,2,3) 


4VY 1 




m2 ) 


2s2 


4VY 


fll 


m2 ) 


2s2 



.In(-I) +in(^^)2 - J-n^ _ Li2(l) -\\r.{^f ) + 0(s) 



2 ' 



m 



12 



L_-LXr,{:^)\^[\-x)sc\{W) 



(III. 14) 



with 



1 -x + z 



2ki ' kji 



(III. 15) 



and Co(l,2,3) Cq(1,2,4) we obtain for the singular part 



r virt., IR div. 



lead, colour 



OC^ / 47lA/^\\„Tn f 1 1 
271''^ I m2 



' s2 s m2 



in(: 



(III. 16) 
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In the sub-leading colour contribution only Cq '' (1,2,4) contains singularities. The contribution from 



Cq ' ^ ( 1 , 2 , 4) is given by 

; virt., IR div. 



_ a.v 1 

subl. colour 271 A'^ 



/j^"(x + x-2z-l)5Co'^(l,2,4). 



(III. 17) 



Using ref. |21] 



(x + f-2z-l)5Co'^(l,2,4) 
rr 



471^^2^1+0)2 f 1. . . 1 



m 



, , , , , -ln((o) + -ln2((o)+2Li2(l-(o)-7in+0(£) 



rr 



4nf/y\+(0^ 1 

I 1 - 0)2 £ 



ln{(o) + {x + x-2z-l)sCo' (1,2,4) 



with 



CO : 



4z 



1 + a/1 



4z 



(III. 1 8) 



(III. 19) 



we obtain 



r virt., IR div. 



subl. colour 



27iAr ^ 



x+x-l 



'l 1+0)2 

Tln(o) 

£ 1 - 0)2 



(III.20) 



for the singular part. There are additional IR singular contributions from the counterterm diagrams. 



They will be discussed in section III.4 . Furthermore we note that in the equations above is the 



result in d dimensions. This follows immediately from the general properties of soft and collinear 
limits in QCD. This is the reason why no 'ys-problem' arises from the IR singularities: Given the 
fact that the IR singularities from the real corrections are also proportional to the Born result in d 
dimensions, finite contributions which would be sensitive to the yj -prescription cancel together with 
the IR singularities, as long as one uses the same 75 -prescription for the virtual and the real corrections. 



III.3 Finite contributions 



Given the divergent contributions shown in eq. ( III. 10 ), eq. ( [II.16| ) and eq. ( III.20 ) we obtain the 
following result for the finite parts: 



' virt., fin. 



2 7, virt., fin. 



lead, colour 



I n^'i- n ;,virt., fin. 
+ Cf «g 



subl. colour 



with 



1 J virt., fin. 



lead, colour 



sj -x(\ -;Cg)+zB[x^ + 5x2^-6x2 + 10a: + 5x^2 - 16xx+ 12x- 10x2 -4 + 3x^] 

(l-X-)2 



-4fix„z 



1-x 



+ x„ 



71 " 



(III.21) 
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+z 



1 



1 +4x(l -x) +x{x-2)x + l{l -xyx— - (1 -x)(l +2x(l -x))- 



, , ;c^x + 8x-7xx-6x2 + 24x-22 

85(2^-; — + 



x2+x-4 



+8g 



l-;c 



(2-x)(l-x)x2 



5 



+ <j6-10;c-6(l-;c)|+z 



(l-;c + z)(2-;c)(l-x) {\-x^z)(\-x 
\-\-x — Az Jc^+Jc-4" 



l-x + 4x + g2(l-T)-4 , _ 
X 1 —x 



So(l,3) 



+ |l6x-z g2 + 3gix+B^5x^jc+llx^-55x + 52xr- 13xx^ + 64(1 -x) - 33x(l -x)^ 



1-x 



4x^ -x2(2 + 3x) +x(4 -X- 4x^) + 2 - 10x+ 14x^ - 4- ^- 

1 — X / x^ — 



1 



■4z 



-32Sz2 



(l-xf 
1-x 



1 -x + 



4(2,3) 



— < — 2x^ 2x(l— Xg)^+z 



-3 + 3gi(2-x)-x + 4x- 



8x 
1-x 



(1 —X 1 —X 

10 + 2jc + lOx^ - x^(4 + 3x) - x(2 +x + 4x^) + 4x^ - 8^— - + 4^— • 



■8(l-x„)4-8— 

* X^ X 



1 (x2x + 8x-7xr-6x2 + 24x-22) 



x2-4z 



(2-x)(l-x + z) 



+ 



16z2 



1-x 



So(2,4) 



+2<{ -x-3(l -Xg)3 + (l-x)x^-x(l+x)-+z 



1 8x, 



X 1 —X 

+(2 + 6x-2x^-x(3-x) + 4j— ^+4(l-x)x4-4x(l+x)-) , ^ ^ 

1 — X x"^ X x^ — 4z 

-8(1 -x)|x- ^1.4(1,2,3) -8(1 -x)|x- ^1.4(1,2,4) 



+2z(l-x)<; lg2 + 7 + 2x-^(l+x) 



1 + 



1-x 



+ 8z 



1 ^ 2Xi, 



-2z<! --gi(l-x)x--+3x+-x2 + 4- 2xr-8z 

2 2 2 1 —X 



_1-X (1-X)2_ 
1-X 



bJ(3,4) 



5C^(1,3,4) 



1 + 2 



1-x 



18xic^ - 14x^ - 4x^x^ + 2x^x + 6x - 9xic - 3x^x + 4x'^ 
1 



+2(l-x)(l + 3x^) 



x2-4z 



^Q(2,3,4) 



(m.22) 



and 



4B ^ 



subl. colour 



2<^2(l-x)+z 



x-2- 



l-2x 
1-x 



ZXs 



(1-X)2 



4(1,3) 



+4z <^ - 5(2x^ + 3xx - 5x - xx^ + 3 - 2x) + ItB (x - x)xg^ - 



2 (l-x)(x-2x) 



{x-2)xg \ , 



1 -x„-4z 



+2^ 2x^-^+z 



4-X + 2 



■4- ^ 



1-X 



2<j2(l+Jc)-x-4^^ ^ 



+4 
_4J _ 



x^-x{x-2){2 + {x-2)-)-2{l-Xg 



1 ^ (x — 2x)(l—x 



x{x^ - 4z) 



X I x^ — 4z 

4(3,4) 



s!,(2,4) 



— 2(1 — 



1 (x-2)xg 



■+z 



-z(l -Jc)|g2 + 1 +4x + 3jc + 2' 



1 -Xg-4z 



(2,4) 



1 -Xg-4z 
x — x 



<(1,3,4) 



+2z| -x(l +x) +4(1 -X,) +4z^ - 



^CS''i(l,2,3) 



+8(;c + ^-2z-l)<^;c- 



2x„z 1 —si 



l-x 



sCo' (1,2,4) 



+2z|2 
+SzB 



1-X 



Xg(l +Xg) — (1 — x)x — 



x(4 — x) — X — 2 



x(l — x) 



1-X 



1-X 
1-Xg 
1-X 



+ 2 



1 -x„-4z 



+2z|2(l-x) 
-SzB 



-2- — 4^) + - + 



(1-X-) 
(x — 2x)(l — x) 
x^-4z 

sC'o{2,3,4) 

1 



1 — x„ 



1 



+ TZ -+4- 



1-x 1-X 2 
l+x„-x- z + 

* 1-X 1-X 



(1— x)^ 1— X x(l— x) (1— Xg)x (1— x)^ Xg_ 

Xg{\-X) 



\—Xn 



l—x ^x—2 ^ 1— X 

4-x-x— ^ + 2- + 2x- -2 

l—x l—x (1— x)^ 



1 - x„ - 4z 



5CS''^(2,3,4) 



-2z 



+ 



+ 



4z 



, r x{2 + x.:}{2-x) ^ l-x 

-6 + 5x+— ^ + 2x- 

l—x 1 — X 

(l-x)(2-3x) 



1-X 

l—Xg 



I —Xff 



x-6 + 4x- 



1-x 



x(4 — x) — X — 2x 



1-X 



1-X 



1-X 



l-x„-4z 



*<''(!, 3,4) 



+2 2(x-xg)(l-xg) 



+z 



1-X 



14(1 - x) - 2x{4 - Xg) - 4x- + Sz-^ 

* 1-X 1-X 1- 



x(l — x)(x — x) 
x„-4z 



+2|2(x-x)(l-x)-2x + z 

x(4-3x) x(3+x) 



2<' 13-8x-2x-Xp^-x-^^ 



1-X 



1-X 



^+4z.,|(2-x-)B + ^||- 



'(1-x) 

1-X 
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l-x 



A-6x + x{A-x)+r -A ^Ix 



XI -xf 



l-x 



l-x 



J l-Xg-Az. 



^ ^d=6,sl,2 



- — { aA{zs)+b,Ml,3)+b23B'oi2,3) +Z.24Bo(2,4) +^34^0(3,4) + C234Q(2,3,4) 



(in.23) 



with g\ defined in eq. ( |III. 12| ), 
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x^-Az 



( 12( 1 - x) - Sjc-' - Gxjf- + I3x'^- 20x + I6xx) , 



(III.24) 



and B defined in eq. ( |[I.12| ). Tlie coefficients a, Z713 , . . . of tlie loop-integrals in the last line of eq. (|IL23|) 



are the same as the coefficients of the corresponding integrals in eq. ([11.22). 



III.4 Renormalization 



In figure we show the Feynman diagrams which arise from Lct{^R,Aj{,m[i,gR) in eq. (III.l). The 






Figure 2: Counterterm diagrams to Z ^ QQg- 



'mass counterterm' shown in the diagrams a, d is given by 

/((Zip — 1)^ — (Zo — l)ms) = /(5Z>i<^ — hZoniRj = i8Z^i{l^ — m«) + /(5Z>i< — 5Zo)nis, (III.25) 

with the momentum counted in the direction of the fermion flow. The first term will result in contri- 
butions proportional to the tree amplitude, while the second term will generate a new structure which 
is not proportional to the tree amplitude. This term should not contain infrared divergences. Inserting 
the explicit results for SZip, and 5Zo it can be easily checked that this term is indeed IR finite. The 
other contribution from the remaining counterterms (diagrams b, c, e, f) are proportional to the tree 
amplitude. Using the relation Zg = {Za)^^^^ |36] valid in the background field gauge as long as one 
uses the MS or MS scheme for the gluon field and coupling, we find that the gluon renormalization 
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enters only through the residue Za of the renormalized gluon propagator. The contribution from the 
counterterms is thus given by 



ct, a+d 



Ct, a+d 



(111.26) 



with T'^'' ^'^^ the contribution from the insertion of /(5Z>i< — 5Zo For the squared amplitude we get 



25Zvp/j^° - 327ia,(A^^ - l)(5Zvi- - 5Zo)5/j6. 



(111.27) 



In the on-shell scheme we have 



5Z°" 



a. 



m 



and 



-rrCf , 
471 \47I/i^ 



a. 



471 



m 



471^2 



12 

- + 4^ + £ 

£ 8iR ' 



-+4K0(£). 



(111.28) 



(111.29) 



To distinguish between infrared and ultraviolet singularities we have introduced £ir by J = 4 + 2£ir. 
This allows us to check separately the cancellation of the UV and IR singularities. To finish the renor- 
mahzation procedure we must include the contribution from the residue of the gluon propagator : 



{zT)^^' 



(111.30) 



which amounts to 



l{{zT)'l^-i)\'r 



tree 1 2 



(111.31) 



at the level of the squared matrix element. (In the on-shell scheme the residue Z'^ of the renormalized 
quark propagator is one by definition.) We have 



47i3r(i-£) 



— (2n'.-ll^)-2Vln^ 



+ 0(£) = l-6z^^ (111.32) 



where «y- is the number of massless flavours and the m, are the masses of the massive quarks, including 
the flavour of eq. (ILl), that contribute to the gluon self energy. We can now write down the entire 
contribution from UV renormalization: 



hi 



UV 



(28Z^" - 5^^ )/26 - 327ia,(A^^ - 1)(5Z^" - 5Z(J")5/26 



m 



hiY^ -2Aa]{N^ -l)CFrY 



m 



471^^ 



£ 



LO 

6 



+ 



(2?!/ - 1 1^) ln( ^ ) + ^ £ In - 8Cf I - 32a^ {N^ -\)CFhhe. 

(111.33) 



471 



Comparing this result with eq. (|IL10|) it follows immediately that all the UV singularities cancel. In 
addition it can be seen that the factor in front of the UV singularities cancels together with the 
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singularities. As we will show in section [V.3| the same is true for the ^-dependence in front of the IR 
singularities. The remaining /^-dependence is thus given by 



a 



1 



^iPoln(£2)C with Po=3(llA^-2n/) 



(III. 34) 



as it should be according to the renormalization group equation. However this is only true as long as 
we renormalize the mass parameter in the on-shell scheme. 

We close this subsection with some remarks concerning the treatment of 75. Since we perform the 
computation of hg, in d dimensions, we have to give a prescription how to treat 75 . We use the following 
substitution in the axial vector current (remember that // is a 4-dimensional index) [32, 33]: 

CX.v 



Y,y5^zr^e,p,8W, with z'^ = I -^Cp + o{a',). 



3! 

The finite 'counterterm' (Zj* — 1) induces a contribution 
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7,LO 
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(in.35) 



(III.36) 



This term is necessary to restore the chiral Ward identities, which can be checked in the present case 
by verifying the following relation: 



hl^{m = 0)=hf{m = 0). 



AVr 



(III. 37) 



Here, the index yA(AV) denotes the contribution to /jg from the interference of the Bom (1 -loop) vector 
current with the 1-loop (Bom) axial vector current. We have checked that after inclusion of the term 
given in eq. (III. 36) the relation shown in eq. (111.37) is indeed satisfied. 



III.5 Checks of the calculation 



The finite parts of the loop contributions to the virtual corrections are given by the rather complicated 
expressions eq. ( [11.22 ) and eq. ( [II.23| ). It is important to perform some checks to ensure their cor- 
rectness. To this end it is useful to realize that the decomposition in terms of (coefficients) x (one-loop 
integrals) is unique in the sense that in the massless case the coefficients of the loop-integrals can 
be uniquely reconstructed from the coefficients of specific logarithms even after having inserted the 
analytic expressions for the loop-integrals. 

Most of the coefficients are highly constrained due to the general properties of infrared finiteness and 
renormalizabihty. As already mentioned above, only the A and Bq integrals contain UV singularities. 
Therefore the fact that all UV singularities cancel after inclusion of the counterterms is an excellent 
check of the coefficients of the one- and two-point integrals, which are each quite involved but combine 
to give the rather simple expression eq. ( III.10| ). As discussed in section [II.2, all IR singularities 
reside in three-point integrals. They are given explicitly in eq. ( [11.16 ) and eq. ( 111.20 ). Therefore the 
coefficients of the IR divergent Co integrals are constrained by the requirement that the IR singularities 
in the virtual corrections cancel against the IR singularities from the real emission processes. We will 
see in the next section that this cancellation really takes place (cf. eq. ( IV.33 ) and eq. ( ^V.34 )). Note 
that only the coefficients of the IR finite Co integrals are not tested by the consideration above. (The 
coefficients of the box-integrals enter the coefficients of the divergent Cq integrals and are thus also 
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tested.) Since all coefficients of tlie loop integrals are computed simultaneously by one code, tiie above 
considerations provide already a quite powerful check. 

We can, however, perform an additional test by comparing our result in the limit m^O with the known 
massless result for h(, given in references [^, 30|, and [31]. In these references the loop-integrals are 



substituted by their analytic expressions in terms of logarithms and dilogarithms. As mentioned above 
we can reconstruct the coefficients of the loop-integrals from the coefficients of specific logarithms and 
dilogarithms. For example in the massless case the dilogarithm appears only in the box-integrals. As 
a consequence the coefficient of the function r{x,y) defined in eq. (4.11) of ref. [|^] must be propor- 
tional to the coefficients of the box-integrals in our result, if we consider the massless case. Note that 
the massless case cannot be derived by simply taking the massless limit, because the massless limit 
of the loop-integrals does not exist in general. However we can still compare the coefficients of the 
loop integrals with the massless result by setting z = in the coefficients and checking which loga- 
rithm/dilogarithm in the results of ref. [ pO| ] are generated from which loop-integrals if one calculates 
the integrals for massless quarks. We do not need to worry about the additional collinear singulari- 
ties which are present in the loop-integrals for massless quarks, as long as we want to check only the 
coefficients of the loop-integrals. 

In this way we can compare our result for example to the one given in eq. (4.9) and eq. (4.10) of ref. 



[gOp. The definition of H(,j in ref. [ |30[ ] differs by a factor of 2 from ours due to a different convention 
in the definition of the electroweak couplings. Taking this into account, we find complete agreement of 
our result with the coefficients of the In^, In and r{x,y) functions in eq. (4.10) of ref. [||]. We also can 
reproduce their rational function, but only up to a multiple of the Born result. This is not surprising, 
since the massless limit induces additional divergent and finite terms and since we use a different 
method to isolate the IR divergences in the real emission processes. For the leading contribution in 
the number of colours we find agreement of our coefficients with those of eq. (4.9) in ref. [ ]30| ] up to 



an overall factor (—1). This is most probably caused by a typographical error which is also present in 



ref. [|29|], but not in ref. [31] with which we fully agree. 



IV Singular contributions from real emission 

In this section we discuss the cancellation of the soft and collinear singularities appearing in the 
loop corrections in particular eq. ( III. 16 ), eq. ( [II.2C ) and eq. ( III.33| ). To cancel them we must in- 



clude the singular contributions from the real emission processes, in particular e^e^ QQqci and 
g+e^ QQgg, with q denoting a light quark. The reactions e^e^ —>■ QQQ'Q' and e^e^ —>■ QQQQ 
do not yield a singular contribution. As a consequence these contributions can be integrated numeri- 
cally in 4 dimensions. In e^e^ QQQ'Q', QQQQ the collinear singularities which would appear for 
massless quarks are regulated by the quark mass and give rise to In(m^) terms. To a certain extend 
these logarithms cancel in the sum of virtual and real contributions. On the other hand in tagged cross 
sections one encounters in general the situation where uncancelled logarithms remain in the final re- 
sult. For example in e^e^ — > QQqq it can be easily seen that such logarithms will be produced when 
the heavy quark pair is produced by a gluon emitted from the light quarks. Due to the fact that this 
contribution is proportional to the couplings of the light quarks to the Z boson it can never be cancelled 
by the virtual contributions. In the presence of several scales of different orders of magnitude these 
logarithms may spoil the convergence of fixed order perturbation theory and indicate a sensitivity to 
long distance phenomena. (Because of these logarithms it is for example not longer possible to take 
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the massless limit!) The uncancelled logarithms may need a treatment by renormalization group tech- 
niques which goes beyond the fixed order perturbation theory. In particular they should be absorbed 
into the appropriate fragmentation functions. These contributions can also be isolated by experimental 
cuts 



To cancel soft and coUinear singularities several algorithms have been developed in the past [[l^, 42,p3]. 
In this work we use the phase space slicing method which was modified to include mass effects in 



3 jets in references [|^, 21 1 and generalized to treat massive quarks in arbitrary jet cross sec- 
tions in ref. [[T^]. Here we restrict ourself to the singular parts from the real emission processes. The 
finite contributions which can be calculated by numerical integration of the existing leading order ma- 
trix elements will be discussed elsewhere. The idea of the phase space slicing is to split the phase 
space of the real emission contributions into regions where all the partons are hard and not coUinear 
(resolved regions) and the remaining regions where two partons are coUinear or a parton is soft (unre- 
solved regions). The separation into resolved and unresolved regions can be performed by introducing 
appropriate Heaviside functions. In the unresolved regions one can make use of the soft and coUinear 
factorization to approximate the matrix elements. Due to this simplifications it is possible to perform 
the integrations analytically over the unresolved regions. 

To fix the notation and to define as clearly as possible resolved and unresolved contributions we outline 
in the next two sections the phase space slicing method applied to the problem at hand. 



IV. 1 Singular contribution from QQqq 

The matrix element for e+e^ — > Q{kQ)Q{kQ)q{kq)q{kq) is singular in the limit kq\\kq. Defining tjj = 
2k j ■ kj we may introduce the following slicing: 

\T {Q{kQ)Q{kQ)q{kq)q{kq))\^ = [e{tqq - S,^,) +e{s^in- tqq)]\'r {Q{kQ)Q{kQ)qikq)q{kq))\\ (IV.l) 

where Sy^i^ denotes an arbitrary mass scale squared which separates between coUinear and resolved 
regions. The resolved contribution 

dOR = &{tQQ-S^in)\'TiQikQ)Q{kQ)q{kq)q{kq))fdR''{kQ,kQ,kq,kq) (IV.2) 

is free of coUinear singularities and can thus be integrated numerically in 4 dimensions. 
Here dR'^ {kQ,kQ,kq,kq) denotes the four particle phase space measure in d dimensions: 

dR\ku...,k^) = l^ (27r)^6(^) {K - g k) fj (^^^ • (IV.3) 

(rA£ is a symmetry factor.) For small Smm the second term in eq. ( |IV.1[ ) can be approximated using the 
factorization of amplitudes in the coUinear limit: 



®{Smm-tqq)\T{Q{kQ)Q{kQ)q{kq)q{kq))\'' 

®{s^n-tqq)\T {Q{kQ)Q{kQ)g{kq + kq))\h]^Pq<i^s{z) (IV.4) 

too 



where Pg^qqiz) is the Altarelli-Parisi kernel [44] (in conventional dimensional regularization) 
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and z denotes the momentum fraction in the collinear limit: 



k^=z{k^ + ki^). (IV.6) 
Using in addition the factorization of the phase-space measure in the collinear limit [16]: 

dR\kQ,kQ,kg,kg) dR\kQ,kQ,P)dR,^Xl. (kg , k^) (IV.7) 

with 



1 (471^2 A £ 
167l2 r(i-£) 



dRconXk^kg) = —^^-^-^dt,^dz[tqqz{l - zT' (IV.8) 



the unresolved contribution to the cross section after integration over the colhnear region is given by 
dOu{QikQ)QikQ)q{k,)q{kg)) = da{Q{kQ)Q{kQ)g{kg))Cqq (IV.9) 

with 

'471^2 



\ f'dz[z{l-z)]-%q^g{z) 



87l2r(l-£) V^mi 

a, 1 /4Vy 1 i_e r2(l-£) 



27Cr(l -£) V ^rain 7 £ 3 - 2£ r(2 - 2£) ' 



(IV. 10) 



In eq. (IV.9) we have identified the sum of momenta + kg with the momentum kg of the outgoing 
'recombined' gluon. In particular the singular contribution to h(, is given by Cqq x h^^. We get the same 
contribution for each massless quark flavour. This results to an overall factor «y with n'j^ the number of 
light flavours. 



IV.2 Singular contribution from QQgg 

The singular behavior of the matrix element for e^e^ Q{kQ)Q{kQ)g{k\)g{k2) is far more compli- 
cated than that discussed in the previous section. While in e^e~ — > QQqq only collinear singularities 
appear we have now both: soft and collinear singularities. In addition only the colour ordered sub- 
amplitudes show a simple factorization in the soft limit. So we may start with the colour decomposition 
of the amplitude which reads 

'r(efe)(2(^2)5(^i)g(^2)) = (7;„r«,),,,^5i + (r„,r,,),,,g52. (ivii) 

Here cq (cq) is the colour index of the (anti-) quark, and ai, 02 are the colour indices of the gluons. 
In our conventions the colour matrices satisfy Tr[rarfo] = bat- The amplitudes 5i, 52 are the so called 
colour ordered sub-amplitudes. The squared amplitude in terms of the sub-amplitudes is given by 

I \'rmQ)QikQ)g{ki)g{k2))\^ = {N^ - 1) (n{\S, |2 + |52|2) - +S2\A . (IV12) 

colour \ / 

Note that the sub-amplitudes are gauge independent. We can thus study the soft and collinear limits 
independently for \Si\^, \S2\'^ and l^i +5212. The contribution |5i +52|2 is essentially QED-like. It is 
therefore free of collinear singularities in the limit that the two gluons become collinear. The amplitude 
^2 can be obtained from Si by exchanging the two gluon momenta. So let us start with the treatment 
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of ^i- The appropriate invariants for the phase space slicing are tgi, tu and tjQ. We start with the 
following identity: 



= 


[0(fei - 


- ■^rnin) + @{s 


■nin - tQl)] [&{h2 - ^min) + @{Smin " ^12)] 












&{tQi - 




■^111111)0(^22 '^min) 


+ 


0(fei - 


^min ) ® ('^min 


- ?i2)0(?2e ~ ■^min) 1>2 collinear 


+ 


@{tQi - 


■^min)® ('^min 


-tn)@{Smm-hQ) 2 soft 


+ 




-tQi)&{Smm 


- h2)®{t2Q - ^min) 1 SOft 


+ 


0(?ei - 


■5min)0(fl2 - 


■5min)0('5min " ?2(2) 2,2 coUincar 


+ 




-tQi)@{tu- 


Smm)&{hQ - ^mm) QJ coUincar 


+ 


®('^min " 


-tQl)&{tu- 


■5'min)0('5'iiiin " ^22) Q'l coUincar and Q,2 collinear 


+ 


®('^min ~ 


-tQl)&{Synm 


- h2)®{i^mm - t2Q) 1 and 2 soft 



(IV. 13) 

The first line in the decomposition above corresponds to the resolved region: both gluons are hard 
and not collinear to each other. Combined with the squared matrix element this contribution can be 
integrated numerically in 4 dimensions. The next three lines describe the single unresolved regions 
which yield a singular contribution after phase space integration. There are also single unresolved 
contributions when one of the two gluons becomes 'collinear' to the heavy quark or anti-quark. These 
contributions do not yield a singular contribution after phase space integration, because of the mass 
of the quarks. In principle they should be of order s^in. In fact for small enough values for s^i^ 
these contributions vanish, because the quark mass serves already as a cutoff. This can be checked by 
including them in the numerical integration procedure. The last two lines describe 'double unresolved' 
contributions. As long as we restrict ourself to three-jet quantities, these terms do not contribute. The 



above slicing of the phase space differs slightly from the one used in our previous work ||20|, ^ and 



corresponds to the one used in ref. [18]. For the leading colour contributions the slicing ([V.13) is more 



convenient than the one used in references [20, 21 ], because the soft factor (see below) is more compact 
and the separation of soft and collinear regions is simpler, facilitating the numerical implementation. 
We have checked numerically that the two slicing methods are equivalent. 

Let us now discuss the contributions from the soft regions, for example the region where momentum 
k\ is 'soft': 

MOu = ®{Smm - tQ\)®{s^m - hlMhQ " ^min) 1^1 1 V/?^ (^Q, , ^1 , ^2) (IV. 14) 

Using the soft factorization of the squared matrix element 

Ar(Ar2_i)|5^|2L^'g2^/(^^,^j,^2)|T(Q(fce)(2(^e)g(^2))P (IV15) 



with 



4/- 4P 4k^ 
f(k- k k )- --^---^-—^ 

hshj lis ^is 



(IV. 16) 



together with the soft factorization of the phase space measure 



dR'{kQ,kQ,kuk2) ^ dR''{kQ,kQ,k2)-dR,oft (IV.17) 
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(the factor 1/2 is the symmetry factor for the two identical gluons) we arrive at 



8dau = ^S{m,0,tQ2)\'riQ{kQ)QikQ)g{k2))\W{kQ,kQ,k2) (IV.18) 
with 

S{mi,mj,tij) = yJ dRsoftf{ki,h,kj)(d{smm - tis)&{smm - t.y), for kj = mj. (1V.19) 
Explicitly 



Sim,0,tQ2) = SFTT^f^l (^]JimAtQ2), (IV.20) 



271 r( 1 V '^min / V '^min 

where 

,2 



J{m,0,tQ2) = @{tQ2-n?) 



£2 2z2\m^) ^2e\m^) 2^^ '^tQ2 



-l^+l-k(2) + i 

2£2 2s 2^ ^ 



(1V.21) 



The case that gluon 2 becomes soft can be treated in the same way. To finish our discussion of the 
singular contribution from let us discuss the contribution from the collinear region. Here the same 
procedure as in the massless case applies. For details we refer to the previous section or to ref. 1 16]. 
So we quote only the final result: 

©(fgl - ■Smin ) (■Smin " ? 12 ) © (^22 ~ ■^min )\Si\^ dR'^ (kQ,kQ,kY 

~ NCgg\'T{Q{kQ)Q{kQ)g{k,))\^dR\kQ,kQ,k,) (IV.22) 



with 



' ^ ' dz[z{l-z)]-'Pgg^,{z) 



2 167l2r(l-£) \Srnin J £ 

a, 1 /Anu^Vl 



and 



r = (^e)''^ I (^e)" 3(i-£)(4-3£)r(i-£)2 

£ £ 2£(3-2£) r(2-2£) 



(IV.24) 



where zq and zq are defined as follows: 



'fig 'eg 

For |52p we may apply the same procedure as for |5i p. The soft singularities appearing in l^i + 
can be treated in a similar way as the singularities in Si . Since the amplitude + ^2 is a QED-like 
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contribution with massive quarks, it does not induce colhnear singularities. The soft singularities can 
be isolated most easily by using the following identity: 

1 = [0(^min - tQX -tQi)+ @{tQi + fgi - Srain)] [®{Smm - tQ2 -^22)+ ©(^22 + tQl " ^min)] • (IV.26) 

If we expand this expression we obtain 

1 = —@{SmiB—tQl—tQi)&{s^ia — tQ2—tQ2) 

+ &{tQi+tQi-S^in)&itQ2+tQ2-Smin) (IV.27) 

The first term describes a double unresolved contribution: both gluons must be soft to satisfy the 
Heaviside functions. This term does not contribute to a three-jet observable. The next two terms 
describe the situation in which gluon 1 or gluon 2 is soft. It is still allowed that both gluons are 
soft, this explains the minus sign in the first line that avoids over-counting. The last term denotes the 
resolved contribution which can be calculated numerically in 4 dimensions. Using the slicing given in 
eq. ( ^V.27 ) we thus obtain for the soft contribution from |5i +52|: 

[0(^min -tQl-tQy)+ ©(^min " ^22 " ^gl)] ( " ^ ) (^^^ " 1 ) I ^1 + -^2 1 ^dR'^ {kQ,kQ,kx M) 

^ -^S{m,m,tQQ)\'T{Q{kQ)Q{kQ)g{kg))\^dR\kQ,kQ,k,) (IV.28) 
where we have used the factorization in the soft region which now takes the form: 

- ^(A^' - 1)1^1 +^21' -^g'sf{kQM,kQ)\T {Q{kQ)Q{kQ)g{ks))\\ (IV.29) 

Note that we have once again identified the momentum of the remaining hard gluon with kg. The 
explicit result for S{m,m,tQQ) is [21] 



where 

,^2 



J{m,m,tQQ) = ^ - i±^ln((o) + |lln((o) - i ln2((o) - 2Li2(l - (o)| + 0{e) (IV.31) 

^^£l— (0 1— CO-'I^S 2 J 

and CO was defined in ( pi. 19| ). Combining all singular contributions from \Si\^, \S2\^, and l^i +52^ we 
obtain the following result for the unresolved contribution: 

NS{m,0,tQg) +NS{m,0,tQg) +NCgg^g - -J-5(m,m,fgg) I da[Q{kQ)Q{kQ)g{kg)) 



271 r(i 



tQQ + 2rrp- 



'mm 

E 



\ x<jA^[(^) J{m,0,tQg)+(^) J{m,0,tQg)-llgg^g 

/ I L \ '^min / \ '^rain / 

^ /(m,m,?gg)| 



da{Q{kQ)Q{kQ)g{kg)) (R^.32) 



To get the singular (unresolved) contribution to h(, one has simply to replace do{Q{kQ)Q{kQ)g{kg)) by 
Note that both quantities do{Q{kQ)Q{kQ)g{kg)) and /jg*-* are evaluated in d dimensions. As we 
will show in the next subsection the same structure is obtained from the virtual corrections. After the 
cancellation of the IR singularities we thus need da{Q{kQ)Q{kQ)g{kg)) or equivalently /ig^ only in 4 
dimensions. 
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IV.3 Infrared and coUinear finiteness 



We close this section by demonstrating that the singular contributions discussed in the previous sub- 
sections indeed cancel the singularities encountered in the virtual corrections. 



Combining eq. ( pI.16| ), eq. ([11.20), and eq. (111.33) we obtain for the IR singularities of the virtual 
corrections the following result: 



r virt., IR div. 



a, 

271 



■N 



i-iln(%)--^ln(% 



1 Jr 

8 



1 1 1+0)2 1 

+--- An{(}))-- 

N el-or 8 



2Cf-^(2«'^-11A^) 
6 ■' 



f,LO 



(IV.33) 



The singular parts of the real emission processes are given by (eq. ( |IV.9| ), eq. ( |IV.32[ )): 

.2\ E 



r real, unres. 
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div 
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r LO 



where we have used the expansions 



J{m,0,tQg) 
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'eg 



(IV.34) 



_ ln( )-ln(— ) 
o ^Qg ^Qg 

+ (iln2(^) + iln2(^) + 17^2 - ^)8 + 0(82), 
Z tn!> Z tn„ J to 



(IV.35) 



with j{m,0,tQg),j{m,m,tQQ) denoting the order 1 contributions in the expansion. They are given by 



+ &itQg-m^ 



+ ln 



1+ )■ 
fQg. 



(IV.36) 
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(IV.37) 
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Note that we have not expanded the factor to show expUcitly the cancellation of this factor as it was 
promised in section III.4 . Comparing eq. ( [V.33 ) and eq. ( [V.34 ) one sees immediately that the singular 
contribution from the real corrections cancel exactly the IR divergences from the virtual corrections. 



V Summary and conlusions 



In this work we have calculated the virtual corrections to the parity-violating functions F3 and Fg of 
the fully differential cross section eq. ([L2) for e^e^ QQg- We have given results for the UV and 
IR singularities as well as for the finite contributions. We have shown explicitly the cancellation of 
the UV singularities after the renormalization has been performed. In addition we have calculated 
also the singular contribution from real corrections using the phase space slicing method. We checked 
that this contribution cancels exactly the IR singularities in the virtual corrections. Together with our 
comparison to the known results for massless quarks the cancellation of the IR and UV singularities 
serves as an important check of our calculation. 

We have given our next-to-leading order results in terms of one function /jg, from which ^3 and F(, can 
be reconstructed using eq. ([L6)-eq. ( II.10| ). The final result for the contribution from three resolved 
partons, ^^^^ can be obtained from the following formula: 



h, 



3 res., NLO 



h 



virt.,fin.^^rest^^ 



real, unres. 
6 



fin. 



(VI) 



where h 



virt., fin. 



is given in eq. ( |III.21| ). The contribution /j?**' is given by 
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(V.2) 



where gi was defined in ( |III.12| ). The finite contribution from unresolved partons is given by 

1 . 
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real, unres. 



fin. 
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(V.3) 
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The functions j{m,0,tQg), and j{m,m,tQQ) are given in eq. ( [V.36 ) and eq. ( IV.37 ), respectively. Note 



that h 



real, unres. 



fin. 



depend on the method used to combine virtual and real corrections. Furthermore we 



can set (i = 4 in /jg*^ in eq. ( V.3) because all singularities are cancelled. 



The results presented in this paper comprise the whole contribution of three 'resolved partons'. Parity- 
violating three-jet observables can be computed using these results together with the contribution from 
four 'resolved partons'. The latter does not contain IR singularities and can, for a given jet algorithm, be 
calculated numerically in a rather straightforward fashion using known results of the matrix elements 



for the 4-parton final states [15 1. In particular one can determine the 3-jet and 4-jet forward-backward 
asymmetries ApB for massive quark jets to order in the QCD coupling, for instance for b quark jets 
at the Z peak. 
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